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Abstract 

The stability of metallic clusters of sodium (Na) in the octahedral 
cages of iVa-doped fullerites NclqCqo and NauC§o is discussed within 
a Thomas-Fermi model. It is shown that the tetrahedral 7Va4-cluster 
in NclqCqq has an electric charge ~ +2.7 (in electron charge units), 
while the body-centered cubic TVag-cluster in NauC§o is almost elec- 
trically neutral. 

The nature of the metallic clusters is currently receiving a great deal of 
interest. m Metallic clusters trapped in zeolites cages provide the opportu- 
nity of studying them by various magnetic resonance techniques, [@]~[f| and, 
in this respect, the nature of the chemical bonding, the distinction between 
metallic and covalent bondings and the degree of ionicity of the clusters 
are of utmost importance. Recently, the existence of iVa-clusters in the 
octahedral cages of iVa-doped Ceo-fullerite has been suggested from X-ray 
analysis. Similar clusters of calcium (Ca) in Ca-doped C^o has previ- 
ously been pointed out.|7j] Tetrahedral-shaped iVa 4 -cmsters seem to build up 
in the octahedral cages of the (face-centered cubic) fcc-Na 6 C 60 ,^ with the 
Na atoms disposed alternately on the corners of a cube of side 2 A. Similarly, 
body-centered cubes of 9 sodium atoms (iVag-cluster) seem to occur in the 
octahedral cages of the fcc-NauCQo,M with the cube side 3.2 A. This is 
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in contrast with the Ceo-fullerite doped with heavier species, such as A^Cqq, 
A = K, Rb, Cs, where the alkali cations push through the /cc-structure 
of the host pristine C 60 and distort it into a fecc-structure, while acquiring 
a conventional ionic character. || The stability and the chemical bonding of 
isolated, small-size alkali-clusters, like Na- and Li-clusters, have extensively 
been studied. 0~[[[2| Recently, the formation of iVa-clusters in Na & C w has 
been challenged, on computational grounds. [13] Here we present a Thomas- 
Fermi model of ./Va-clusters in No,qC§q and NanCeo, discuss its stability and 
compute the distribution of electric charges and the degree of ionization of 
the clusters. 

The main observation we start from is that the Na-Na distance in these 
./Va-clusters is much shorter than the Na-Na distance 3.7 A in bulk metal- 
lic sodium, though slightly larger than twice the iVa + -radius 1.15 A. In- 
deed, the Na-Na distance in the tetrahedral iVa 4 -cluster in A^a 6 C 6O 0] is 
2^/2 A = 2.82 A, while in the cubic iVag-cluster in ./Van Ceo || this dis- 
tance (the shortest) is 3.2^/3/2 A = 2.75 A. This indicates that, in com- 
parison with the bulk metallic sodium, the iVa-clusters in doped Cqq are 
squeezed (or compressed) in the octahedral cages provided by the neigh- 
bouring fullerene molecules. In the octahedral coordination of the fcc- 
Cqq the fullerene molecules are placed at ~ 7.1 A from the coordination 
centre, [PU'RIB"! and, allowing for a radius of ~ 3.5 A of the fullerene molecule, 
one obtains an average radius R = 2 A of the octahedral cage.||'|| For the 
tetrahedral iVa,4-cluster, whose atoms are placed at R = 2^3/2 A = 1.73 A 
from the coordination centre, this cage provides an average atomic radius 
of 2 s A, and an average inter-electron separation much smaller than the 
Bohr radius. This suggests that a Thomas-Fermi model is suitable for de- 
scribing the iVa4-cluster. The confinement of the cluster in the octahedral 
cage is ensured by the electronic clouds of (carbon) C-atoms on the sur- 
face of the fullerene molecules, which, in order to be penetrated, require 
higher-energetical electrons. They act as very high potential barriers. The 
confinement of electrons in the much narrower tetrahedral sites of the fcc- 
Cqo (whose radius is estimated at ~ 1.12 A J/J) may certainly provide them 
with enough energy as to make them able of penetrating the confinement, 
which may explain the absence of clusters in these cages. Similar considera- 
tions apply for heavier-alkali atoms in the octahedral cages, and even for the 
iVag-cluster in iVaiiC 60 , whose atoms, placed at R = 3.2\/3/2 A = 2.75 A, 
clearly trespassed the 2 A limit found above. However, the penetration of the 
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electronic clouds of the C-atoms proceeds gradually, and, in the octahedral 
cages, the C-atoms provide another limit of about 4 A in the (lll)-directions: 
it is the distance from the coordination centre to the octahedral faces passing 
through the centres of the fullerene molecules. We may take, tentatively, a 
mean value of these two limits (2 A and 4 A), and set up an average cage 
radius Rq ~ 3 A for the iVag-cluster. This gives an average atomic radius 
3s A, which again justifies a Thomas-Fermi model for this iVag-cluster. It is 
worth-remarking at this point that the resistance the cage walls oppose to the 
penetrating electrons (which may also apply in the case of the zeolites cages) 
is an expression of the Pauli exclusion principle, which requires fast ly- varying 
(in space) wavefunctions for the latter (and therefore higher energy) if they 
are going to penetrate into the space regions already occupied by atomic-core 
electrons. This confining effect of the cage boundaries on an atomic cluster 
may be formalized as an external field "pseudo-potential", much in the same 



manner as it has already been done for bulk metals. [|I^j It can easily be seen 
that these "pseudo-potentials" would vary strongly over short distances (as 
in the short-range parts of the inter-atomic Lennard- Jones or Born-Mayer 
potentials, for example), and their functional dependence on distance (which 
would depend, among others, on the number of confined electrons) will pro- 
vide an effective cage radius R for a given cluster, as we shall discuss below 
in more detail. 

As it is well-known the Thomas-Fermi model is a qusi-classical theory 
of interacting electrons and atomic nuclei which assumes that the Fermi 
wavector kp of the electrons, and the electron density n, 

- = ^4, (i) 

vary slowly in space. This is true (except for distances very close to the nuclei) 
due to the effective screening of the nuclear charges accomplished by the high- 
density electrons. We shall see that this assumption is a consistent one also 
in the present Thomas-Fermi model of clusters. Under this assumption the 
particular shape of the atomic cluster is no longer very relevant, so that we 
may introduce an additional simplification by assuming that the total charge 
z of the nuclei is uniformly distributed over a spherical shell of radius R, plus 
a central ./Va-nucleus of charge +11 (in electron charge units) in the case of 
the iVag-cluster. We denote by V the electrostatic potential created by this 
distribution of positive charges, and remark that its particular expression 
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is not very important in subsequent calculations; the only requirement on 
the positive charge distribution is that of symmetry, from obvious reasons of 
stability, and the approximation made is that of radial symmetry. In the high- 
density limit the electron exchange interaction is a higher-order contribution 
(as well as other contributions to the electron interaction energy), and we 
may assume that an electrostatic potential tp is generated by the electron 
distribution, which does not depend on hp. We may therefore set up the basic 
Thomas- Fermi relationship (in atomic units Bohr radius an = h 2 /me 2 = 
0.53 A, e 2 /a H = 27.2 eV and h = 1) 

\k 2 F -<p-V-U = -<p , (2) 

where U is the potential well of the cage and (p is the chemical potential; the 
latter should be constant in order to ensure the local equilibrium and equal 
to the "pseudo-potential" of the external field at the cage frontier Rq, such as 
to prevent any flux of electrons pouring in and out of the cage. The variation 
of this external "pseudo-potential" over the distance R may be taken as the 
depth U of the potential well. Using (1) and (2) we may rewrite the Poisson 
equation Aip = 4irn as 

A(v? + V + U - v? ) = ^(y? + V + U - <p )% , (3) 

everywhere except for r = R; introducing the reduced variables x = r/R 
and 

* +v+u -" = {smfl (4) 

we arrive at the Thomas-Fermi equation 



x 
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x 2 X =X 2 , (5) 

for any x between and 1 except for x = a = R/R . Since n is a continuous 
function <p and its two first derivatives must be continuous functions; it 
follows that the function x has a slope discontinuity (but is itself continuous) 
exactly as the derivative of V has, i.e. at x = a; the magnitude of this 
slope jump will be obtained below by Gauss law. It is therefore convenient 
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to define Xi( x ) = x( x ) for < K a and Xi{ x ) = x( x ) for a < x < 1, with 
the continuity condition 

Xi (a) = X2 {a) • (6) 

The number of electrons inside a sphere of radius x is easily obatined 
from (1) — (5) as 



N(x) = AuRl /* dx -x 2 n = ^ /* dx 



x X -X 



9tt 2 



^Xi-Xi + Xi(0) , 0<x<a, 



(7) 



9tH 



i28m 



f ) { x X2-X2 + a[x' 1 (a) -x 2 (a)} + Xi(0)} , a < x < 1 , 



and one remarks that it is continuous at x — a. The total number of electrons 
in the cage is obtained from (7) as 



N 
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{x' 2 (l) - X 2 (l) + a [ X [(a) - x' 2 (a)} + Xi(0)} 



(8) 



By Gauss' law we get easily the net charge inside the sphere of radius x as 
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q{x) = 



128R 3 



(x~ x x) , 



(9) 



whence, by using (7), we obtain the positive charge 

9tt 2 

qi+ix) = 128i?f l(0) 

inside the sphere of radius x < a. We remark that, if it exists, it is concen- 
trated at x — 0. Obviously, for the iVa4-cluster we should require 



(10) 



Xi(0)=0 



:n) 



while for the body-centered cubic iVa 9 -cluster the corresponding boundary 
condition is 
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128.Rg 



X i(0) = 11 . 



(12) 



The total charge inside the sphere of radius x, a < x < 1, given by (9) 
consists of z-positive charges distributed over the spherical shell of the cluster, 
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a charge +11 at the centre, if the case, and —N(x) electron charges; whence, 
by using (7), we get the magnitude 
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128.Rg 



XiO) - xi(°) 



(13) 



in the slope discontinuity of the function Xi as expected; (13) provides the 
third boundary condition for our problem, with z = 44 for the iVc^-cluster 
and z = 88 for the iVag-cluster. The fourth boundary condition needed for 
the uniqueness of the solution of (5) is provided by the stability condition of 
the positive shell of the cluster. It is easy to get the electric self-field acting 
on the shell; in the absence of the central atom it is given by 



E s = z/2R 2 , 
while in the case of the centered cluster we obtain 



(14) 



E s = z/2R 2 + 11/ 'R 2 . (15) 
The electric field on the shell due to the electrons is readily obtained as 



E r = - 



dip 



9vr 2 1 



Xi(a) - ax'i{a) - xi(0) 



(16) 



dr 1 a 128^ a 2 

(or equivalently in terms of X2), and we remark that, in contrast to (14) and 
(15), it is negative, i.e. the electrons try to stabilize the exploding cluster 
shell. Using (13) we arrive easily to the equilibrium condition 



2xM 



(17) 



both for the non-centered and for the centered cluster. The equation (5) is 
solved numerically with the boundary conditions given by (6), (11) — (13) 
and (17), both for the Na^- and the iVag-cluster, and the results are described 
below. In each case we have been interested especially in the total charge 
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128i$ 



Xa(l)-xi(l) 



(18) 



in the cage, as given by (9). 

In the case of the iVa4-cluster we have used R = 1.73 A and Rq = 
2 A (a = R/Rq = 0.86), as suggested by the X-ray diffraction data. 
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The function x is obtained by solving numerically the above equations; it 
corresponds to a total charge q = +2.72, i.e. an average charge +0.68 per 
each iVa-ion. As we have discussed above, the numerical results obtained 
for the function x are valid everywhere except for a range of about twice 
the Bohr radius around the position of the positive shell of charges, where 
the variation of x/ x is too large; even so, one can estimate that almost 
half of the electrons are located there, ensuring thus the screening of the 
positive charges. Another limitation of the present calculations occurs from 
the assumption of radial symmetry for the function x- We may estimate 
the error introduced by this approximation as follows. Assuming the same 
variation of the electron density per unit length in all directions the relative 
contribution of the angular part of the variation of the potential derivative 
with respect to the radial one is 5ip' a /5'(p r ~ \f^, f° r a cluster consisting of 
s atoms; on the other hand, as these small variations are proportional to the 

variations of the distance we have (<Vr) + 2 = i^'Pr)^ wnere the 

latter corresponds to the neglect of the angular part. As the variations of the 
potential derivatives are also proportional to the charge (from the Poisson 
equation), we find then easily that the relative error made in estimating 



the charge is at most ~ 1 — yl + 8n/s/(l + 2yAn/s). In the case of the 
iVa4-cluster this amounts to~ 40%; one can see that the radial-symmetry 
approximation improves its accuracy for large s. 

Increasing the amount of negative charge on the cluster the position of 
the maximum value of x moves toward smaller values of x (smaller values 
of R), and it exhibits an opposite behaviour on increasing the degree of 
ionization. The variation of x with respect to the cage radius R brings into 
discussion the energy of the cluster. One can check easily that the density of 
the electronic kinetic energy is given by 



Ekm ~2^k dk - k '- ( 10j ^io ) - \ ' ( 19 ) 

where (2) and (4) have been used; whence one can obtain straightforwardly 
the total kinetic energy of the electrons 



^-^1^4- (20) 

£2 



20R 7 Jo 

We obtain similarly the potential energy of the electrons 
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£ po t = -^0- fdx-K- N<p ■ (21) 
so that the electronic energy may be expressed as £ e \ = £\ — N<p , where 

e^JMff ix .4. (22) 

The self-energy of the positive charges is given by £ s = £ — zU for the 
non-centered cluster and £ s = £q — (z + 11) U for the centered one, where 
£ = z 2 /2R in the former case and £ = z 2 /2R+ llz/R in the latter. On the 
other hand, expressing the interaction energy of the shell with the electrons 
in two distinct ways, namely, zip (R) = — J dr nV, we get 

U -Vo = ^x'Al) • (23) 

Summing up all these contributions we obtain the total energy of the cluster 

as 



£ = £ Q + £ 1 + £ 2 -{2z-q)^ , (24) 

where £ 2 = —z (U — ip ) as given above; in the case of the centered cluster we 
should replace z by z + 11 in (24). The numerical integration gives E\ = —28 
and £ 2 = —42 for the A^a 4 -cluster, while £ = 296; the total energy, consisting 
of (24) and the energy of the electrons transferred to the cage walls, must 
be negative, in order the cluster be formed. For small values of the total 
charge of the cluster the latter contribution to the total energy is small, so 
that we may estimate that the energy E given by (24) must acquire negative 
values. Hence it follows that tpo must be of the order of (296 — 28 — 42) /N ~ 
2.6 (for N = 44). On the other hand, ipo equals the external "pseudo- 
potential" at the cage frontier R , as we have discussed before, whence one 
may have an estimation of the magnitude of this "pseudo-potential". One 
can see, therefore, that it acquires high values, as to ensure the confinement 
of the squeezed cluster within the Rq range. In addition, as the present 
estimations are valid over a scale length larger than the Bohr radius, one can 
also see that the external "pseudo-potential" has a high variation over short 
distances. This illustrates again the gradual resistance developed by the 
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inner-core electrons of the environment to the penetrating "pressure" of the 
confined cluster. The electron density varies also fastly over short distances at 
the frontier of the cage, in order to ensure the delimitation of the compressed 
cluster from its confining environment. The "pressure" exerted by the cluster 
electrons on their environment results in the deformation of the electronic 
clouds of the surrounding atoms. Having determined (fo by the external 
"pseudo-potential" at the cage frontier R , the two remaining parameters 
i?oand R will be obtained, within a complete theory, by minimizing the total 
energy of the cluster (given by (24)) and the electrons transferred to the cage 
walls. Of course, this requires first the computation of the "pseudo-potential" . 

From the discussion above one can see that the results of the computation 
are very sensitive to the values of R and Ro; this sensitivity may result in 
a higher inaccuracy for clusters consisting of a greater number of atoms. 
This can be seen in the case of the iVag-cluster, where, according to the 
experimental data suggested by the X-ray analysis,^] we have used R = 
2.75 A and Ro = 3.15 A (a — R/Rq — 0.87). The function \ is obtained 
again by solving the corresponding equations given above; it corresponds, 
within errors of ±0.03 A in R and Rq, to a total charge of the cluster q ~ ±1. 
At this level of approximation we may conclude that the centered-cubic Nag- 
cluster is practically neutral. The same discussion regarding the screening 
and the value of the chemical potential (fo (~ 2.4) applies here as for the 
-/Va/i-cluster. It is worth-noting that a similar calculation for a non-centered 
cubic cluster consisting of 8 atoms of Na, gives a total charge q ~ +1, for 
the same values of R and Ro. Finally, noticing that the theory contains the 
charge z only in the product zRq (see, for example, (12), (13) and (18)) 
we can say that similar clusters formed by lithium may also form in the 
octahedral cages of the fullerite with a smaller size, or, equivalently, with a 
diminished cluster-size and ionization charge. 

In conclusion one may say that the calculations presented here indicate a 
metallic bonding of the iVa 4 - and iVan-clusters squeezed into the octahedral 
cages of the iVa-fullerides. We emphasize that this squeezing is a central 
feature of the picture presented here, whose direct consequence is that of 
affecting all the electron states in the cluster; enough electronic energy may 
be gained in this way as to stabilize the cluster at lower values of R, as 
compared with R . This is in contrast, at least for the iVa4-cluster, with 
recent claims[0 that the iVa-atoms in Na§C§o are fully ionized and that the 
octahedral iVa-ions are placed much further away from the coordination cen- 
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tre. These claims, based on numerical simulations within the local- density 
approximation, make use of a certain distinction between the valence and 
the core electrons, which might be inappropriate for the squeezed clusters. 
As we have seen from the present calculations, almost half of the electrons 
seem to be affected by the presence of the neighbouring atoms in the clus- 
ter, which indicates a drastic change in the inner electronic states of the 
isolated atoms. The metallic character of the alkali clusters in fullerides may 
be tested by various techniques, most notably by Raman and 23 iVa-NMR 
spectroscopies. [|17]] In addition, it is worth- mentioning that the metallic clus- 
ters squeezed in atomic cages may exhibit two types of collective modes, 
associated, one, with the vibrations of the positively-charged ions and, an- 
other, with the oscillations of the electron density. Preliminary calculations 
indicate that these two modes are in the range of the ultra-violet optical 
spectroscopy and the electron spectroscopy. 
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